
tic
%% Pre script
close all

%% Config
Is  = 18.8e-21;
n   = 2;
V_d = 1:.02:2;
V_t = 0.0258;
inc = .1;

%% Calculate

i=(exp(V_d/(n*V_t))-1);


%% Plot
figure
hold on;
for step=0:10
    y=Is*ones(1,size(V_d,2));
    plot3(V_d,(step*inc)*y,(1+step*inc)*Is*i);
end
xlabel('Volts');
ylabel('I_s');
view(-45,45);
grid on
hold off;
toc